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Abstract 

We use character sums to confirm several recent conjectures of 
V. I. Arnold on the uniformity of distribution properties of a certain 
dynamical system in a finite field. On the other hand, we show that 
some conjectures are wrong. We also analyze several other conjectures 
of V. I. Arnold related to the orbit length of similar dynamical systems 
in residue rings and outline possible ways to prove them. We also show 
that some of them require further tuning. 

1 Introduction 

In a recent series of papers, V. I. Arnold [H IS |3l HI [5], [6] has considered dy- 
namical systems related to linear transformations in finite fields and residue 
rings and made a number of conjectures. We observe that the study of the 
length, distribution of element and other properties, of the orbits of such dy- 
namical systems has a long and successful history, which dates back to early 
works of N. M. Korobov [46], H. Niederreiter [511 [60], A. G. Postnikov ^ 
and many other researchers. Here we show that some classical results im- 
mediately imply some of these conjectures. We also show that several other 
conjectures are not correct as they are stated in [H El El HI El E] and need 
some adjustments. 
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For a prime p and a positive integer n, we denote by Wpn the finite field of 
p"^ elements (we refer to [S2] for the background information on finite fields). 

We fix a primitive root of Wpn and recall that we also have Wpn = Wp [d] . 
In particular, Wpn can be considered as an n-dimensional vector space over 
Wp, where with each element a G Wpn one can associate the coordinate vector 
a = (oq, . . . , a„_i) G Wp from the expansion 

n-1 

j=o 

Accordingly, V. 1. Arnold [1] suggests to study the sequence of vectors a^ = 
(ao,m, • • • , a„_i_m) £ Wp corresponding to the powers 

n-1 

^- = J2<',,-^'-'- (1) 

j=0 

Clearly, assuming that Wp is represented by the elements of the set 
{0, 1, . . . ,p — 1}, one can view the points 

-a^, m = 1, . . . , M, (2) 

p 

as M points of an n-dimensional unit cube [0, 1]". For M = p" — 1 these 
points form a regular cubic lattice (with only one missing point (0, . . . , 0)). 
It has also been conjectured by V. I. Arnold [1] that in fact even the first 
M < — 1 powers already form a rather uniformly distributed point set. 
Namely, given a region Q G [0, 1]" with smooth boundary, we denote by 
N^{M,Q) the number of points ([2]) which belong to Q. The conjecture of 
Section 2. A of [1] asserts that 

N^{M, n) = Mvol n + o{M) (3) 

provided that M ~ /ip" for some fixed /i > (and p oo). 

We start with an observation that using classical bounds of incomplete 
exponential sums with exponential functions, see [lH |l6l [52] , and some stan- 
dard tools from the theory of uniform distribution, see Section 12.41 one can 
derive the following improved version of the conjecture ([3]): 

iV^(M, n) = Mvol Q + (Mi-i/>i/2nQQgp)i+i/n^ ^ 
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which is nontrivial whenever M/p^/^(logp)"'"'"^ ^ oo. In fact, using some 
results of H. Niederreiter [SSI ED] one can easily extend the above result in 
several directions. 

In fact, using the results of J. Bourgain and M.-C. Chang [13j, which in 
turn generalize several recently emerged results of J. Bourgain, A. A. Gli- 
bichuk and S. V. Konyagin one can also study the distribution in 

intervals of the set for extremely small values of M. For example, see [TB] 
for more details and a version of the bound (jl]) which is nontrivial provided 
that M > for any fixed e > and sufficiently large p. The bound of the 
error term in [16] is not completely explicit, so for large values of M the 
bound (jl]) is better than that of [TBj . 

Moreover, for n = 1, that is, for prime fields, using bounds of exponential 
sums from [9l [SH |36], one can obtain nontrivial results for even smaller 
intervals, which however holds only for almost all primes p (rather than for 
all p) . 

Furthermore, motivated by the results of [11], |32l, we consider the dis- 
tribution of vectors ci^ where instead of an initial segment [1,M], m runs 
through the values of a polynomial. Unfortunately, we are not able to treat 
arbitrary polynomials with integer coefficients for every primitive root d but 
rather obtain a result which holds for almost all primitive roots. However, 
in the case of monomials, employing the bound of exponential sums with the 
sequence , m = 1,2, . . ., from [27|, we obtain a nontrivial estimate for 
every primitive root 

V. I. Arnold [H [S] [S] |6] also describes similar dynamical systems in the 
residue ring modulo £ and makes several conjectures about the length 
of the orbits. More specifically, given an integer g > 2 with gcd{g, €) = 1, 
V. I. Arnold [U E] El E] suggests to consider the dynamical properties of the 
residues g"^ (mod i). 

We recall that the Carmichael function X{i) is defined for all £ > 1 as the 
largest order of any element in the multiplicative group Z^. More explicitly, 
for any prime power p", one has 

\(n-\ = l ifp>3orz/<2, 
^yP ) <^ 2^-2 if p = 2 and 1/ > 3, 

and for an arbitrary integer i > 2, 

A(£)=lcm(A(pr),---,A(prO), 
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where i = . . .p''/ is the prime factorization of L Clearly, A(l) = 1. We 
also let ip{€) denote the Euler function, which is defined as usual by 

s 

^(£)=#z*=nprfe-i)' 

with V5(l) = 1. Finally, for a an integer g with gcd((7,£) = 1, we denote by 
tg{C.) the multiplicative order of g modulo C.. Clearly, we have the divisibilities 

t,(£) I m I 

Several conjectures of [H [3l O [6] can be reformulated as various statements 
about the relative size of tg(£), A(£) and ^{d), on average and individually. 
We discuss these conjectures and show that some of them are already known 
in the literature, while some can be proved to be wrong. It is suggested in 
Section 1 of [5] that ioi g = 2 the average multiplicative order 



grows like 



L 

1=1 

gcd(g,£)=l 



T,(L)~cM— (5) 

for some constant c{g) > depending only on g (we note that in [5] it is 
made explicit only for g = 2). 

We show that the classical result of Hooley [12] on Artin's conjecture, 
implies, under the Extended Riemann Hypothesis , that the conjecture ([5]) is 
wrong and in fact 

Tg{L) > i^exp (C((7)(logloglogL)3/2) . 

for some constant C{g) > depending only on g. Furthermore, we believe 
that in fact Tg{L) grows even faster. It is possible that the method of proof of 
Theorems 1 and 2 in [8J, which in turn is an extension of the method of [SZ] 
(see also together with the result of Hooley can be used to derive 
that, under the Extended Riemann Hypothesis, 

> p^exp ((logloglogL)2+°(i)) . (6) 

log L 
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For the upper bound it is probably natural to assume that 

Note the sum on the right hand side of ([7]) has been estimated by P. Erdos, 
C. Pomerance and E. Schmutz 

Finally, we give a guide to the literature concerning results and methods 
which can probably be of great use for the theory of algebraic dynamical 
systems over finite fields and rings. 

It is very well known that there are close ties between number-theory and 
dynamical systems. For example, one can associate dynamical systems with 
continued fractions, various number systems, the 3x + 1 transformation and 
other number-theoretic constructions. A wealth of very interesting results 
can be found in the literature. 

However, we would like to use this paper as an opportunity to attract 
more attention of the dynamical system community to a great variety of al- 
ready existing number theoretic results and techniques which can be of great 
significance for studying various algebraic dynamical systems. In particular, 
these include, but are not limited too, bounds on various exponential sums, 
periods of various sequences and average values of associated arithmetic func- 
tions. For this very purpose we do not try to formulate and prove our results 
in their full generality but rather limit ourselves to the most interesting and 
illuminating special cases. We however indicate possible extensions of our 
results and directions for further research. 

Although exponential sums have been used for this purpose, see the work 
of M. Degli Esposti and S. Isola [18] and of P. Kurlberg and Z. Rudnick ^8] . 
their full potential seems to be not fully used in the dynamical system theory. 
We would like to stress that all such applications follow the same pattern: 

bounds of exponential sums 

ij- 

distribution in aligned boxes 

^ 

distribution in arbitrary regions with smooth boundary 

^ 

ergodic properties of the corresponding dynamical system 
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The link between exponential sums and the distribution in aligned boxes 
is provided by the Koksma-Sziisz inequality, see Theorem 1.21 of [TU] . 

The link between the distribution in aligned boxes and arbitrary regions 
is given by the results of H. Niederreiter and J. M. Wills [64j and their more 
recent refinement of M. Laczkovich |39] . 

Surprisingly enough, the essentially tautological link between the distri- 
bution in arbitrary regions and ergodic properties has never been exploited in 
a systematic way, although it definitely deserves much more attention which 
we hope to attract with this paper. 

2 Dynamical systems in finite fields 
2.1 Preliminaries 

Here we show how well known bounds of exponential sums can be used to 
derive various results about the orbits of 'd-^^"^\ with a polynomial f{X) e 



Throughout this section, any implied constants in the symbols O may 
depend on n and Q (and occasionally, where obvious, on an integer parameter 



As we have mentioned, the results of in this section can be extended in 
several directions. 

2.2 Background on finite fields 

Let Uq, . . . , ujn-i be a basis of F^n over Wp which is dual to the basis 1, "(9, . . . , ■(?' 
That is. 



k). 




0<i,j<n — 1 



where 



n-1 



a' 



is the trace of a G Wpn in Wp. 
Therefore, from ([1]) we derive 



Tr 



(8) 
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for every j = 0, . . . n — 1 and m = 1, 2, . . .. 

It is also useful to recall that there are ip{p'^ — 1) primitive roots of F^n. 



2.3 Background on exponential sums 

Let us denote Gp{z) = exp{2'7iiz/p). Then for every 7 G F^n the function 
a I— >■ (7Tr (a)) is an additive character of Wpn. 

For example, it follows immediately from a combination of Theorem 8.24 
and Theorem 8.81 of [52] (see also [45^ i46j and the references therein), that 
for any M < — 2 the following bound holds 



max 



M 

E 

m=l 



ep(Tr(7^™)) 



0(p"/2logp). 



(9) 



The following estimate is a special case of a more general result of 
We remark that in [35] it is shown only in the case n = 1 but the proof 
extends to arbitrary fields without any changes. 

Lemma 2.1. For any primitive root eW pTi J any two subsets X, y G 2pn_]^ 
and any function ip{y) with 



max 

y&y 



the following bound holds 



max 



^V^(y)ep(Tr(7r^)) 



yey 



O i'^p 



,9n/8+o(l) 



3/4 



We know recall the bound of exponential sums with from [27J which 
we use in the proof of Theorem 12. 7[ More precisely, we use Theorem 6 (for 
k = 2) and Theorem 7 (for A; > 3) of [27] (we remark that in [27] these results 
are proven only for n = 1 but the general case can be obtained by a simple 
typographical change of p to p"'). 

Let us define 



1 



p{k) = { 



\k/2] - 1 
L 2k \k/2] +2' 



if = 2; 



if > 3. 



(10) 



7 



Lemma 2.2. For any primitive root G IFpn, the following bound holds 



max 



m=l 



< p(l-p(fc))n+o(l)_ 



2.4 Background on discrepancies 

For a finite set W C [0, 1]" and domain Q C [0, 1]", we define the fl- 
discrepancy 

#{uGWnfi} 
A{U,n)= — — --vo\n 

and the box discrepancy of U, 



D{U)= sup A(W,B), 

BC[0,1]" 

where the supremum is taken over all boxes B = [ai, x . . . x [ak,(3k]- 
We define the distance between a vector u G [0, 1]" and a set F C [0, !]"■ 

by 

distfu, F) = inf llu — w|| 

wGr 

where ||v|| denotes the Euclidean norm of v G IR". Given £ > and a domain 
VL C [0, 1]" we define the sets 



fi+ = {u G [0, I dist(u, Q) < e) 



and 



fi; = {u G ^] I dist(u, [0, < . 

Let b{e) be any increasing function defined for £ > and such that 
lim£^o^(^) = 0. Following 091 [61], we define the class of domains 
n C [0, l]'^ for which 

vom+ < 6(e) and yoIVL^ < b{e). 

As special case of a result of H. Weyl [69] implies that for domains Q with 
a piecewise smooth boundary, one can take b{e) = 0{e). 

Lemma 2.3. For any domain Q C [0, 1]" with a piecewise smooth boundary 



voinf = 0{e) 
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A relation between D{IA) and A(Z^, Vt) for f2 G M.h is given by the follow- 
ing inequality from [49] (see also [M]). 

Lemma 2.4. For any domain Vt G M.h, we have 

A{U,n) =0{h {n^/^D{Uf'^)) . 

The Koksma-Sziisz inequality, see Theorem 1.21 of [T^], provides an im- 
portant link between box discrepancy and exponential sums: 

Lemma 2.5. For any integer L > 1, and a set U C [0, l]" of M points, one 
has 



( 



D{U) = O 



where 



n-1 



}_ J_ T-r 1_ 

M ^ J-J- (1 + |c,-|) 

:(co,...,c„_i)GZ"\{0} j=0 ^ I 
\cj\<L j=0,...,n—l 



J]exp (2 



TTtC ■ U 



n-1 



CjUj 



j=0 

denotes the inner product of c = (cq, . . . , c„_i) and u = (mq, . . . , Wn-i) ■ 

To estimate the box discrepancy of the set ([2]) we apply Lemma 12.51 with 
L = [p — l)/2. By ([8]) we see that the corresponding exponential sums takes 
shape 



M 



'n-1 



M 



j=0 



5^e,(Tr(7^-)) 



m=l 



m=l 



where 7 = Cqujq + . . . + Cn-ioJn-i ^ lF*n. 

Applying the bound (Q together with Lemma 12.51 we see that the box 
discrepancy of the set ([2]) is O (M~^p"'/^(logp)""'"^) , see also [SHI EO] and 
references therein for several more general results. Now the bound (jl]) follows 
directly from Lemmas 12.31 and 12.41 



2.5 Distribution of points in orbits 
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For an polynomial f{X) G "^[X] and a given region Q G [0, 1]" with 
smooth boundary, we denote by N^{f; M, Q) the number of points 



P 



m 



(11) 



which belong to Q. 



Theorem 2.6. Let f{X) G 2[X] be a fixed nonconstant polynomial and let 
be the set of all ip{p"' — 1) primitive roots ofWpn. For any positive integer 
M < p^ — 1 and any region Q G [0, 1]*^ with piecewise smooth boundary, we 
have 

^ |Ar^(/; M, fi) - Mvol n\ < M^-^/ 



1?G0 



Proof. To estimate the box discrepancy of the set ffTTj) we apply Lemma 12.51 
with L = (p — l)/2. As in Section [2^ by ([H]) we see that the corresponding 
exponential sums takes shape 



M 



'n-1 



M 



j:e,[Y: c,Tr ^^/(-)) ]=J2e, (Tr (7^/^)) 

j=0 



m=l 



m=l 



where 7 = cqujo + . . . + Cn-iUJn-i G 1F*„. Applying Lemma [2751 we see that 
the box discrepancy D^{f; M) of the set ([2]) satisfies 



Y^D,{f-M) = 0[- + 



{logpY 



max 

p M 7GlF;n 



M 



5^e,(Tr(7^9^(-))) 



M] 



Let X be the set of all elements of Zpn_i which are relatively prime to — 1. 
Fix an arbitrary primitive root G F^n. Then = {i^g | x G X}. 

We denote by 3^ the value set 3^ = {f{m) (mod — 1) | m = 1, 
and by ip{y) the multiplicity of ?/ G 3^ (that is, the number of m = 1, . . . , M 
with y = f{m) (mod — 1)). In particular, < M and by the famous 
Nagell-Ore theorem (see [13] for its strongest known form) we have 

\1/ = max |V^(y)| = p°^^\ 

y&y 



We derive from Lemma 12.11 that 

M 



max 



^e,(Tr(7^^('"))) 



m=l 



^ ^9n/8+o(l)j\^3/4^ 
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which imphes the bound 



Mf-, M) < /"/8+«(i)m-i/^ (12) 



Now, since Q has a piecewise smooth boundary, from Lemmas 12.31 and 12.41 
and the Holder inequality, we derive 

J2 W&if-^ M, n) - Mvo\n\ = oImJ2 Di){f] M)i/" J 
= O I M i^{p- - l)"-i g M) j 



O (m^(p"- 1)1-1/" M)j 



l/n^ 



Since k/(p{k) = 0{\og\ogk) for every integer k, see Theorem 328 of 
from fll2l) we derive the desired estimate. □ 

For example, we see from Theorem 12.61 that for every fixed e and M > 
pn/2+e^ for almost all primitive roots of '(9 G Fpn, we have N^{f;M,iT) ~ 
Mvoin for every region Q G [0, 1]" with smooth boundary. 

For primes p such that — 1 has a certain prescribed arithmetic structure 
M.-C. Chang [16] obtained nontrivial results which hold for all primitive 
roots, rather than on average. Such primes are rather sparse but one can 
show that there are infinitely many of them. 

Now, for an integer k > 1 and a given region Q G [0, 1]" with smooth 
boundary, we denote by N^^k{^) the number of points 

-a^fc, m = 1, . . . - 1, (13) 
P 

which belong to Q. 

As before, we define p{k) by ( ITOl) . 

Theorem 2.7. For any primitive root i} G Wpn and any region Q G [0, 1]" 
with piecewise smooth boundary, we have 

N^,ki^) = p"vol n + (p«-'^(^)+"(i)) . 
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Proof. Arguing as in the proof of Theorem I2.6[ and using Lemma 12.21 instead 
of Lemma I2.H we easily deduce that the box discrepancy of the set ( TT^ 
can be estimated as O . Now applying Lemmas 12.31 and 12.41 we 

conclude the proof. □ 

We remark that similar results can be obtained in a more general situation 
(for example without the request that ^ is a. primitive root). We however 
follow the settings which exactly correspond to those of Section 2. A of Jl]. 

2.6 Some other conjectures and open questions 

In Section 2.B of [1] a conjecture is made which essentially means that the 
consecutive values 'd"^, -(9™-+^ are independently distributed. It is easy to 
see that this is incorrect. For example, if n = 1 and -(9 = 2 is primitive 
root modulo p (see [12]) then if ■(9™' G (p/4,p/2) then ^9™ G {p/2,p), which 
happens for about p/4 values of m = 1, ... ,p — 1, while a conjecture given 
in Section 2.B of |5 predicts that this should happen for about 3p/16 values 
of m. 

Sections 2.C and 2.D of [4J contain a number of interesting questions 
about the geometric properties of the set of points ([2]). We remark that 
the bound O (M~^p"'/^(logp)"^^) on the box discrepancy of ([2]) obtained in 
Section 12. 4^ immediately implies that any aligned cube [a, inside of the 
unit cube [0, 1]" with the side length (3 — a > CM~^/"p^/^(logp)^+-^/", for an 
appropriate constant C > 0, contains at least one point ([2]). This immedi- 
ately implies upper bounds of the same order on the largest distance between 
the points ([2]) and on the largest radius of a ball inside of [a, which does 
not contain any points ([2]). Moreover, using some standard modifications, 
see [17], one can drop the logarithmic factor from these bounds. 

3 Dynamical systems in residue rings 
3.1 Preliminaries 

Given an integer g > 2 with gcd{g,i), V. I. Arnold [H El El E] suggests to 
consider the dynamical properties of iterations of the map x i-^ gx in the 
residue ring (which is equivalent to studying the residues g"^ (mod i), in 
particular to studying the multiplicative order tg{i)). In particular, in the 
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papers PQ El El E] a number of suggestions have been made about the average 
orbit length of this and several similar dynamical systems. 

We remark that indeed if the orbit length is sufficiently large then, fol- 
lowing the standard scheme, one can derive some analogues of (jl]) from well 
known bounds of exponential sums |l5l HHl E2l ES 160] . 

Here we provide a brief guide to the literature and demonstrate that many 
existing techniques are suitable for studying these questions and in fact imply 
that some conjectures of [H El El E], based on numerical calculations, need 
some further adjustments. 

Throughout this section, the implied constants in the Landau symbol 'O' 
and in the Vinogradov symbols '<^' and may occasionally, where obvious, 
depend on and are absolute otherwise (we recall that U <^V and V ^ U 
are both equivalent to the inequality U = OiV)). 

3.2 Analytic number theory background 

Let T^gix) denotes the number of primes p < x, such that g is a. primitive 
root modulo p. 

We recall the following celebrated result of Hooley [12] : 

Lemma 3.1. Under the Extended Riemann Hypothesis, for every integer g 
which is not a perfect square, there exists a constant A{g) > such that 

X 

7ig{x) ~ A{g)- . 

logx 

Let tt{x\ k,a) denote the number of primes p < x with p = a (mod k). 
We need the following relaxed version of the Brun-Titchmarsh theorem, see 
Theorem 3.7 in Chapter 3 of [39]. 

Lemma 3.2. For any integers k,a > 1 with 1 < k < x, the bound 
holds. 

Let V be the set of prime numbers. The following estimate can be de- 
rived via partial summation from Lemma [3.21 see, for example, the proof of 
Theorem 3.4 in 1221 . 
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Lemma 3.3. For any integer / > 1 the bound 



E 



1 log log X 



•C 



^-0,2^^ plog(3x/p) ip{f)\ogx 

p€V, f <p<x 
p=l (mod /) 

holds. 

Proof. Let h = \ 2 log f\,H= [log x] . Then 

1 ^ 1 

P.vl^<p<. ^'l^g^Sx/p) ^^^^^^^^^^.^^ plog(3x/p) 

p=l (mod /) P=l (mod /) 

« log(3xe-^-i) ' ^{f)j ^ ^{f) ^ j (log(3x) - j - 1) ' 

The result now follows. □ 



3.3 Average multiplicative order 

Here we show that the conjecture ([5]) is wrong and in fact Tg{L) grows faster. 

Theorem 3.4. Under the Extended Riemann Hypothesis , for every integer 
g which is not a perfect square, there exists a constant C{g) > such that 

Tg{L) > r^exp (C(^7)(logloglogL)3/2) . 
log^ 

Proof. Let Vg be the set oi p eV for which g is a, primitive root modulo p. 

Let us put 

Q = exp ( A/log L). 

For an integer k > 2 we consider the set Cg{k,L) of positive integers i G 
[L/2, L] of the form i = pi . . . p^ where pi G Vg and Pi > Q, i = 1, . . . , k. 
By Lemma |3.H considering only those integers i G Cg{k,L) for which 
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Pi, ■ ■ ■ ,Pk-i < {L/2QY/^^ ^\ and thus L/2pi . . .pk-i > Q-, we have 



#C,{k,L)>- 



E 



E 



> 



{A{g) + 


[m 


2k\ 




iA{g) + 




2k\ log 


L 


{A{g) + o 




2k\ 



Pl,...,Pk-l€'Pq Pk&Vq 

Q<pi,...,Pfc_i<(L/2Q)i/(''-i) -£'/2pi...Pfe_i<Pfc<i/Pi 

1 



E 



pi,...,Pk^ie'Pq 

Q<pi,...,Pfc_i<(L/2Q)i/(fc-i) 



Pi . . .pk-i\og{3L/pi ...pk- 



E 



pi,...,Pk-ieVq 

Q<pi,-,Pfc-i<(W2Q)i/"--^) 



Pi ■ --Pk-i 



( 



\ 



k-l 



E - 

PGP ^ 
VQ<p<{L/2Q)i/(fe-l) / 

By partial summation, we derive from Lemma 13.11 that 



V - ~ A((7)(loglog(L/2g)i/(^-i) - loglogg) ~ ^loglogL, 

^ P 2 

Q<p<{L/2Q)i/('=-i) 

uniformly for k with log A; = o(loglogL). Therefore, uniformly for k with 
logfc = o(loglogL) 



#£,(fc,L)> 



(A(^) + o(l))^L(loglogL)^ 
2'=A;!logL 



We now denote by Qd{k, L) the set of positive integers i < L of the form 
- Pi . . .pk where Pi eV are distinct primes, > Q, i = 1, . . . , fc, and 



max gcd(pj — l,pj — 1) = d. 
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For d < Q^l"^ we use Lemma 13.21 to derive 



{k-2)\ 



EE El 

pi,...,Pfc_2eP pev gGV 
pi,--;Pk-2>Q P>Q Q<g<L/ppi...pk-2 

p=l (mod d) q=l (mod d) 



Pl,...,Pk-2&'P 
Pl,-,Pk-2>Q 



E 



1 



pev 

Q<P<L/pi--Pk-2 
p=l (mod d) 



ppi... pk-2 \og{3L/ppi . . . pk-2d) ' 



Applying Lemma 13131 with f = d (once) and theh / = 1 (/c — 2 times), we 
obtain, that for d < Q}!"^ 



#Q,(A;,L)< 



L(log logL) 



fc-i 



{k-1)\Lp{df\ogL' 



Since ^{d) ^ (i/ log log (i, see Theorem 328 of [40j, we conclude that for any 
L> > 0, 

,L(loglogL)''^^ 



Q^/^>d>D 

Also using the trivial bound 



-l+o(l); 



(A; -2)! log L 



Pl,P2e'P, 
Pl,P2<L 
pi=P2=l (mod d) 



P1P2 



E - 



n<L 



^2 



^71=1 (mod d) j 



for d> Q^/'^ we derive that 



E md{KL)<Q-"^Li\ogL)' 

d>Q^I'^ 



Hence, for D < Q^'"^ 



Y*Qci{k,L)<D- 



_i+o(i)^(loglogL) 



fc-i 



d>D 



(fc-2)!logL 
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Finally, let S{L) be the set of positive integers i < L such that for some 
p > Q. Clearly 

#5(L) < J2 [L/P'\ « L/Q. 
per 
P>Q 

Therefore, for Dk — ?,^, for the set 

\d>Dk 

we have 

^ 2^fc!logL ' 

provided log/c = o(loglogL). 

On the other hand, for every i — pi . . .pk G T^giL) we have 

tg{e) = \Cm{tg{pi),...,tg{pk))=lCm{pi-l,...,Pk-l) 

k i—1 ^ 

s nfe-i)ng,dta-i.ft-i) 



Therefore, under the above condition on A;, we derive 



^ #7^,(A•.I) ^ + o(l))^-I(loglogI)' 



ee'Rg{k,L) ^ 
gcd{g/)=l 

L(loglogL)*^exp(0(A;3)) 
logL 

Taking 

/c= c(c/)yiogloglogL 

for an appropriate constant c{g) > 0, depending only on g (which guarantees 
that the term (log log L)'^ exceeds, say, the square of factor exp(0(A;^))) we 
finish the proof. □ 
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As we have remarked, we believe that the bound of Theorem 13.41 is not 
tight and in fact a stronger bound can be derived by using the method 
of [57], modified in a similar way as that of [8] to deal only with special 
primes, see also p5] . 

It is a very interesting question to obtain more precise information about 
the behaviour of Tg{L), for example to establish whether ([7]) is correct. It 
is possible that the method of [23] combined with the methods and results 
of [3Zl[5T] are able to handle this task. In fact, even already existing results 
of [211 HH [51], without any modifications or adjustments, may shed light on 
many issues risen by V. I. Arnold in [H [3l [5l [6] . 

We also remark that a dual question about the the average value 

g=l 
gcd(g,£)=l 

is studied in [STJ [531 EB]- particular, in [56] one can also find various 
upper and lower bounds on T(£), while its behaviour on special sequences is 
considered in [371 113] • 

It is well known that if an integer g > 1 is fixed then for any function e{x) 
with e{x) — *■ as x — *■ cxo, for almost all primes p the bound tg{p) > p^/^+^^p) 
holds, see [23l [26l [65] for various improvements of this result. For almost 
all integers i, similar type bounds are given in [47j . 

It is clear that when g varies, tg{i) runs through divisors of A(^). In 
fact through all the divisors of A(£). Accordingly the question about the 
behaviour of r(A(£)) becomes of interest (where T{k) is the number of all 
integer positive divisors of k > 1). Partially motivated by this relations, 
F. Luca and C. Pomerance [55j obtained tight bounds on the average value 
ofr(A(£)). 

3.4 Average additive order 

V. I. Arnold [5J also asks about the average period Q{i) of the map x ^— x + a 
in taken over all a G Z^. This function, which can be expressed as the 
following sum over the divisors of i 

d\e 
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has been studied in detail in [3Z1 IMj . 

In particular, it is shown in Theorem 3.1 of that 

^E^(^) = ^^ + 0((logL)^/='(loglogLr/^) 

where C,{s) is the Riemann (-function. This gives a more precise and explicit 
form of the assertion made in |5] that on average Q{i) grows linearly. Clearly, 
for any prime p, we have T{p) = Q{p — 1), thus it is also interesting to study 
Q{i) on this and some special sequences of £, see [37]. More results on 
arithmetic properties of Q{i) have been established by F. Luca |54j . 



3.5 Average divisor 

We note that well known bounds of number theoretic functions implies that 
the assertion made in [5] that the "average divisor" of ip{i) is 



Tim) logc 

(where cr(/c) is the sum of all integer positive divisors of A; > 1) with some 
f3 = 0.96 ± 0.02 is false. As it follows from the classical number theoretic 
bounds 

2 < T{k) < 2(i+°(i))i°g'=/i°si°s'=, k + l< a{k) = O(fcloglogfc), 

see Theorems 317 and 323 of [IQ], respectively, for any sufficiently large k we 
have 

. ^ . ....... 

Finally, we mention that the suggestion made in [5] that the average value 
of d{k) behaves like 



' ' 2\ogK 

k=i ^ 

is wrong too. Since d{k) is a multiplicative function, so is d{k)/k, which also 
satisfies the conditions of the Wirsing theorem, see [70]. Thus one can easily 
show that in fact 
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for some absolute constant k = 0.4067. . ., see (TU] for this and some other 
results on the properties of the average divisor, including an asymptotic 
expansion of D{K). However, the question on the average value of d{(p{i)) 
is harder and is of ultimate interest. 

It can also be relevant to mention that average divisor d{k) takes integer 
values for almost all integers k > 1 but is almost never a divisor of k, see pUj 
ESIEB] and the references therein. Certainly similar questions about d{(p{i)) 
are very natural and interesting. 

4 Repeated squaring and other nonlinear 
transformations 

Using [71 |12l |271 |28l [291 EI] one can also easily derive various uniformity of 
distribution results for the vectors where e > 2 is a fixed integer. Alter- 
natively, these results can be interpreted as results about orbits of repeated 
powering x i— > x'^. In particular, with e = 2 one can study the distribution 
of elements in orbits of repeated squaring x i— > in finite fields and rings, 
see [21 [3] where the corresponding dynamical system is outlined. The results 
of [271 [2H1 [291 [31] show that if the orbit is long enough then the vectors 
a^m are uniformly distributed. Using the bound of J. Bourgain [12], one can 
consider very short orbits (and not necessary fixed values of e) although the 
bounds obtained within this approach are less explicit. 

These results are complemented by the estimates on the orbit lengths 
of such transformations which are obtained in |30l [58] and which show that 
these orbit lengths tend to be large (and close to their largest possible values). 

The distributional properties of dynamical systems generated by general 
non-linear transformations x i— » f{x) where / is a rational function, over a 
finite field or a residue ring, have been extensively studied in the literature 
as well, see [20l [211 [HI [621 [HSj and the references therein. As in the case of 
repeated powering all these results indicate that if the orbit is long enough 
then its elements are uniformly distributed. On the other hand, these results 
are still missing their essential counterpart, namely estimates on the orbit 
length. Obtaining such estimates (for general or specific functions /) is a 
very important open question. 
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5 Further remarks and extensions 



Results of a different flavour but also describing the distribution of powers 
of primitive elements modulo a prime p are given in |67j . 

We have already remarked that analogues of our results hold for an ar- 
bitrary not necessarily a primitive root, provided the multiplicative or- 
der of is large enough. For more general formulations of Lemma 12.11 and 
Lemma [2^2} see [331 135] and [27], respectively. 

Analogues of the bound ([9]) and Lemmas 12.11 and 12.21 are also known for 
residue rings 2^, see [l5l HG], [52] and [281 129], respectively. Thus one can 
study orbits of i)^ in residue rings as well. In particular, it is well known 
that one can use these bounds to obtain various uniformity of distribution 
results suggested in Section 2 of [3]. 

Furthermore, the dynamical system corresponding to the repeated squar- 
ing of a unimodular matrix has been considered in [2|. Using the results and 
methods of [5B], one can prove various uniformity of distribution properties 
of orbits of such dynamical systems. Accordingly, the results of [3 SH \5U\ can 
be used to obtain similar statements for analogues of the above dynamical 
systems on elliptic curves over finite fields. 
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